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Abstract. We describe an attempt to model numerically Bose-Einstein correlations (BEC) from 
"within", i.e., by using them as the most fundamental ingredient of some Monte Carlo event 
| generator (MC) rather than considering them as a kind of (more or less important, depending on 

O j. the actual situation) "afterburner", which inevitably changes original physical content of the MC 

D ' code used to model multiparticle production process. 
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Introduction. The problem of BEC is so well known that we shall skip introductory 
remarks (referring in that matter to other presentations at this conference and to [1] for 
(Nj . more information) and instead we shall start right away with our subject concerning 

the proper numerical modelling of BEC, which we call BEC from within. Although 
phenomenon of BEC is with us from the very beginning of the systematic investigation 
of multiparticle production processes, its modelling is still virtually nonexisting. With 
the exception of attempts presented in [2, 3] all other approaches are tacitly assuming 
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^3 ' that on the whole BEC constitute only a small effect and it is therefore justify to add it 
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in some way to the already known outputs of the MC event generators widely used to 



model results of high energy collisions (in the form of the so called afterburner) 1 . There 
J|h \ are two types of such afterburners: 

(a) those modifying accordingly energy-momenta of identical secondaries (and cor- 
recting afterwards the whole sample for energy-momentum conservation) - they 
apply to each event separately; 

(b) those selecting events which already have (due to some fluctuation present in any 
MC code) right energy-momenta of identical secondaries and counting them (by 
introducing some weights) many times (see [1] for references) - in this case energy- 
momentum balance is left intact but instead the particle spectra provided by MC 
code are distorted (albeit in most cases only slightly); they apply only to all events. 

Example of modifications introduced by afterburner. It must be stressed that modifi- 
cations that such afterburners introduce to physical backgroundbehind a given MC code 
were never investigated. Tacit assumption made is that they are small and therefore irrel- 



Methods presented in [2] were never used in practice, only [3] applied their approach to e + e annihila- 
tion data. 
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evant 2 . The problem is that with the MC codes in use at present it is practically impos- 
sible to estimate the nature and strength of such changes. In can be done only by using 
some simple scheme of cascading, for example simple cascade model for hadronization 
developed by us some time ago [4]. This model per se leads to no BEC (cf. left panel 
of Fig. 1). To obtain effect of BEC one has to add to it some simple afterburner, for 
example the one proposed by us in [5]. In it one preserves the whole space-time and 
energy-momentum structure of each event but changes the charge assignment to sec- 
ondaries in such way that clusters of identically charged particles occur (the original 
number of (+), (— ) and neutrals remains the same). This automatically leads to BEC 
3 . The right panel of Fig. 1 shows what are the changes introduced by this afterburner 
in the original cascade: the multicharged vertices occur now (with total charge in the 
whole event remaining conserved). It means that, in principle, one could obtain BEC 
directly from MC (i.e., without any afterburner) allowing in it for the appearance of 
such vertices (according to some prescribed scheme - for example as some multicharged 
clusters, such possibility has been already mentioned in [6])). However, it would be then 
extremely difficult to run such cascade till the very end without producing spurious mul- 
ticharged particles not observed in nature. Our afterburner can be then considered as a 
kind of shortcut realization of such possibility with well defined physical consequences. 
We argue therefore that each afterburner changes original MC code it is attached to in a 
similar fashion but this statement cannot be at the moment substantiated 4 . 

BEC from "within". The above observation was one of our motivations to look for 
the MC scheme, which would be build around BEC rather than starting from some 
single-particle observables. This idea has been for the first time used in [3] where par- 
ticles were selected from a grand-canonical-like distribution with temperature T and 
chemical potential /i chosen in such way as to describe rapidity and multiparticle distri- 
butions 5 and particles were then distributed into some rapidity cells of given width, each 
cell containing only particles of the same charge; the method was then (quite success- 
fully) used to describe annihilation data. Introduction of these cells was the real 
origin of the BEC observed. Actually the idea that BEC demands that particles of the 
same charge are emitted from some cells (named elementary emitting cells - EEC) was 



2 Actually the frequent practice is to use the original MC data to obtain single particle spectra and to 
calculate the corresponding BEC afterwards by means of some afterburner. 

3 Notice that such reassignment of charges results in the effective bunching of momenta very much similar 
to that assumed in the method (a) mentioned above [1]. 

4 Notice that the problem, which is clearly visible in the CAS model, is not at all straightforward in other 
approaches. However, at least in the string-type models of hadronization, one can imagine that it could 
proceed through the formation of charged (instead of neutral) color dipoles, i.e., by allowing formation of 
multi(like)charged systems of opposite signs out of vacuum when breaking the string. Because only a tiny 
fraction of such processes seems to be enough in getting BEC in the case of CAS model, it would probably 
be quite acceptable modification. It is worth to mention at this point that there is also another possibility 
in such models, namely when strings are nearby in the phase space one can imagine that production of 
given charge with one string enhances emission of the same charge from the string nearby - in this case 
one would have a kind of stimulated emission discussed already in [7]. 

5 Actually, in [3], which was using information theory approach based on Shannon entropy, T and jx are, 
in fact, two lagrange multipliers obtained from energy conservation and charge conservation constraints. 
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FIGURE 1. Example of charge flows in MC code using simple cascade model for hadronization [4]: 
left panel - no effect of BEC observed; right panel - after applying afterburner described in [5] (based on 
new assignement of charges to the produced particles) one has BEC present at the cost of appearance of 
multicharged vertices. 



proposed earlier in [8] and consist also cornerstone of our present proposal, which can 
be regarded as generalization of that presented in [3]. This is nothing else but attempt to 
numerically realize the bunching of particles as quantum statistical effect used already 
in connection with BEC long time ago [7] and is done by dividing all available energy 
among particles (taken here as being pions) in such way that a number n ce u of EEC's, 
each containing particles of the same charge, is formed, with multiplicity of particles in 
each EEC followinig Bose-Einstein (or geometrical) distribution 6 . When energy distri- 
bution from which particles are selected is thermal-like then P{n ce u) is Poissonian and 
the total multiplicity distribution is of Polya-Aeppli type [9], closely resembling Nega- 
tive Binomial distribution obtained in the so called clan model [10] (which differs only 
by the fact that particles in clans are distributed according to logarithmic distribution, 
not geometrical one [1]). Particles in a given EEC can have energies spread around the 
energy E\ of the first particle defining this EEC with some width o. With such energy 
spreading allowed one gets quite reasonable results for C2(Qinv) distributions (see [1] 
for details) 7 . 

Here we would like to present extension to this algorithm, which in addition to 
bunching accounts also for the symmetrization of the two-particle wave function (not 
used before) and allows to obtain in addition to CiiQinv) aiS0 CiiQxyz), i- e -> m a sense it 
is 3-dimensional extension of our algorithm. This extension is based on the observation 
that symmetrization correlates the energy-momenta of particles with their space-time 
locations. The bunching of particles considered before was done only in the energy- 
momentum space and left us with a number of EEC's, each with a number of particles 
with well defined energies, Ej (and momenta pi = (Ef — m 2 ) 1 / 2 ). Each EEC is build 



To this end to the first particle selected in a given EEC one adds (up to first failure after which new EEC 
is selected) another particles with probability P = Pq- e~ E l J ', where Po € (0, 1) is constant, E is energy 
of the first particle and T parameter (corresponding to temperature in thermal models). Such form of P 
ensures the characteristic Bose-Einstein form of energy distribution. 

7 Actually, as was shown in [11] (see also [12]), this spreading is crucial to obtain the proper shape of C2 
function. 



Bose-Einstein correlations from "within" 



February 2, 2008 



3 



2.0 

1.9- 

1.8- 

1-7- 

1.6- 

1-5 

1.4- 

1.3- 

1.2- 
1.1- 
l.O- 

0.9- 



T=3.5 GeV, c = o.3*T GeV, P=o.7»exp(...) 

W= 46.5 GeV 

- - W = 91.2 GeV 
W = 182.4 GeV 



1-3 



0.2 0.4 0.6 o.S 1.0 1.2 1.4 1.6 1.8 2.0 
[GeV] 



<y 1.1 

w 

o 



0.9 





•■.1. 




. !'■ 

-si. 

■..-;'>^ 


lfl V : ':-. 




ikTA.,..- 









-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



1 



2.0 
1.9 
1.8 

1-7 
1.6 
1-5 
1.4 
1-3 
1.2 
1.1 
1.0 



0.9 



• a o = 0.1 

a o = 0.3 > T - 3.5 GeV, P=o.7*exp(..J 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

a [GeV] 




FIGURE 2. Example of results obtained with 3 -dimensional version of our model. Upper panels show 
results for different mass of the hadronizing source. Lower panel demonstrates dependence on the allowed 
energy spread in EEC. 



around some "seed" particle, which is taken as particle i = 1. This was enough to get 
Ci{Qinv) (see [1]), but not C2(<2 W ) involving components of p t , Pi( x ^ z y To assign 
them one has first choose some space-time positions for particles in a given EEC taking 
them from some distribution function p(r,t). Actually in what follows we shall use 
only static source approximation, i.e., hadronization is instantaneously and therefore 
p(rj) = p (r) . Now Pi( x ,y,z) have to be correlated with the corresponding space positions, 
H = (xi,yi,Zi), in the way emerging from the symmetrization of the wave functions 
resulting (for the plane wave approximation) in the famous 1 + cos(<5p • 8r) expression. 
Technically this is achieved by accepting only such momenta Pi( x ,y,z) > which for given 
r, lead to cos(<5p • 8r) < 2 ■ Rand — 1 where Rand is random number uniformly chosen 
from interval (0, 1). 

In Fig. 2 we present examples of our new results (the elder results can be found in [1]) 
obtained with full, 3-dimensional version of our model for p(r) being sphere of radius 
R = 1 fm and assuming that all Pi( Xj y jZ ) are spherically symmetric. As one can see now 
in addition to C2(Qinv) one has also corresponding to it Pi( X y yZ ) C2(Q x . y . z ). 

So far we are assuming direct pion production. However, the inclusion of Coulomb 
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and other final state interactions in our approach is straightforward - one must simply 
change cos(- • •) term arising from plane wave approximation used to form obtained by 
some distorted wave function. One can also easily include resonances and allow for finite 
life time of the emitting source. Finally, so far, for simplicity reason, only two particle 
symmetrization effects have been accounted for. Namely, in a given EEC all particles 
are symmetrized with the particle number 1 being its seed, they are not symmetrized 
between themselves. This seems to be justified because majority of our EEC's contain 
only 1 — 2 or 3 particles. But to fully account for multiparticle effects one should simply 
add other terms in addition to the cos(. . .) used above. This, however, would result in 
dramatic increase of the calculational time 8 . Nevertheless the effect of including at least 
terms when symmetrization between, say, particles 2 and 3 are added to the already 
present symmetrization between 1 and 2 and 1 and 3, must be carefully investigated 
before any final conclusion is to be reached. 

Summary. To summarize, we are proposing numerical scheme of modelling quan- 
tum statistical phenomenon represented by BEC occurring in all hadronization pro- 
cesses. Distinctive features of our scheme not present in other propositions are: 

• identical particles are emitted from EEC's and only these particles are subjected to 
BEC; 

• inside each EEC particles are distributed according to the geometrical (or Bose- 
Einstein) distribution; 

• altogether they show characteristic Bose-Einstein form of distribution of energy. 

As result we obtain a kind of quantum clan model with Negative-Binomial-like mul- 
tiplicity distributions and characteristic shape of C2{Qi m ) function [1] (notice that we 
automatically include in this way BEC to all orders given by the maximal number of 
particles in a given EEC). To get also Ci{Qx,y,z) one has to use some additional space- 
time information and the characteristic 1 +cos(<5r • bp) correlations between space-time 
and energy-momenta induced by the symmetrization of the respective wave functions. 
So far this is only a case study, we cannot yet offer any attempt to compare it with ex- 
perimental data. On the other hand our approach offers new understanding of the way in 
which BEC are entering hadronization process. 

We shall close with remark that there are attempts in the literature to model numer- 
ically BE condensation [13] (or to use notion of BE condensation in other branches of 
science as well [14]) using ideas of bunching of some quantities in the respective phase 
spaces. 



In fact, this converges to the proposition presented long time ago in [2], the only hope is that in our 
case symmetrization is performed within given EEC and therefore the number of terms of cos(. . .) type 
involved is rather limited, whereas in [2] the whole source had to be symmetrized at once resulting with 
number of terms growing like n\ where n is observed multiplicity. 
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